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Abstract. A calculation of the current-quark-mass-dependence of nucleon 
static electromagnetic properties is necessary in order to use observational 

^ I data as a means to place constraints on the variation of Nature's fundamental 

parameters. A Poincare covariant Faddeev equation, which describes baryons 
as composites of confined-quarks and -nonpointlike-diquarks, is used to calcu- 

^ ', late this dependence. The results indicate that, like observables dependent on 

the nucleons' magnetic moments, quantities sensitive to their magnetic and 
charge radii, such as the energy levels and transition frequencies in Hydrogen 
and Deuterium, might also provide a tool with which to place limits on the 
'nI" ' allowed variation in Nature's constants. 

o 

oo 
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1 Introduction 



It is a feature anticipated of models for the unification of all interactions that 
5^ I the so-called fundamental "constants" actually exhibit spatial and temporal vari- 

ation. In consequence there is an expanding search for this variation via as- 
tronomical, geochemical and laboratory measurements [H [2] . An interpretation 
of these measurements can materially benefit from calculations of the current- 
quark-mass-dependence of observables characterising hadronic and nuclear sys- 
tems. 

One example is found in the behaviour of hadron masses jSj S] . A variation 
in light-meson masses modifies the internucleon potential, and a variation in 
the nucleon's mass affects the kinetic energy term in the nuclear Hamiltonian. 
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2 Variation of nucleon magnetic moments and radii 

Such changes modify nuclear binding energies and can thereby have a substantial 
impact on Big Bang Nucleosynthesis [5]. These results enable the use of obser- 
vational data to place constraints on the variation of Nature's constants; e.g., 
Ref. [6]. 

Other effects of a variation in hadron masses are an alteration in the location 
of energy levels and in the positions of compound resonances in heavy nuclei. 
For example, Refs. [Tj [8] explore the sensitivity to changes in the light-quark 
masses, m, of the nuclear clock transition between the ground- and first-excited 
states in ^^^Th and the position of the 0.1 eV compound resonance in ^^'^Sm. 
It is noteworthy that the shift of the Sm resonance, as determined from the 
Oklo natural nuclear reactor, currently provides the best terrestrial limit on 
the temporal variation of Nature's fundamental parameters; namely, \Xq/Xq\ < 
1.6 X 10-18 y-i,Xg:=m/ylQCD [ailOllII]. 

Hadron magnetic moments can also depend upon current-quark mass. Calcu- 
lations of this dependence are necessary for the interpretation of measurements 
of quasar absorption spectra and superprecise atomic clocks in terms of the 
variation of Nature's fundamental parameters [12^ \T3[ 114] . It is notable that in 
rudimentary constituent-quark models for hadron bound-states composed of de- 
generate light-quarks, in which the constituents are weakly bound and hence the 
bound-state's mass is accurately approximated as the sum of constituent-quark 
masses, the bound-states' magnetic moments are independent of quark mass. 
Any quark-mass variation in hadron magnetic moments is therefore a gauge of 
hadron structure. 

Herein we report a calculation of the mass-dependence of neutron and proton 
magnetic moments based on a Poincare covariant Faddeev equation model for 
the nucleon. Since the magnetic and electric form factors are obtained simul- 
taneously, we also describe the variation of the nucleons' magnetic and charge 
radii. The background material for our calculation is provided in Sects. [2] and [3l 
and in appendices. Our results are presented and discussed in Sect.Hl Section [5] 
is a brief summation. 



2 Nucleon Model 

In quantum field theory a nucleon appears as a pole in a six-point quark Green 
function. The pole's residue is proportional to the nucleon's Faddeev amplitude, 
which is obtained from a Poincare covariant Faddeev equation that adds-up all 
possible quantum field theoretical exchanges and interactions that can take place 
between three dressed-quarks. 

A tractable truncation of the Faddeev equation is based |15j on the obser- 
vation that an interaction which describes mesons also generates diquark corre- 
lations in the colour-3 channel [16]. The dominant correlations for ground state 
octet and decuplet baryons are scalar (O"*") and axial- vector (1''') diquarks be- 
cause, for example: the associated mass-scales are smaller than the baryons' 
masses [171 [IB] , namely (in GeV) 

^lud]^+ = 0-7 - 0.8 , m(„„)^^ = m(„rf)^_^ = rrndd)^+ = 0.9 - 1.0 ; (1) 
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Figure 1. Poincare covariant Faddeev equation, Eq. HA.11[) . employed herein to calculate nu- 
cleon properties. ^ in Eq. (|A.l|l is the Faddeev amplitude for a nucleon of total momentum 
P = Pq + Pd- It expresses the relative momentum correlation between the dressed-quark and 
-diquarks within the nucleon. The shaded region demarcates the kernel of the Faddeev equa- 
tion, Sect. lA.2l in which; the single line denotes the dressed-quark propagator, Sect. lA.2!T] F 
is the diquark Bethe-Salpeter-like amplitude, Sect. lA.2!2l and the double Ime is the diquark 
propagator, Sect. lA.2!3l 



and the electromagnetic size of these correlations [19] 
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is less than that of the proton. (The last result is an estimate based on the 
p-meson/vr- meson radius-ratio p0t[2T].) 

The kernel of the Faddeev equation is completed by specifying that the quarks 
are dressed, with two of the three dressed-quarks correlated always as a colour-3 
diquark. As illustrated in Fig[Tl binding is then effected by the iterated exchange 
of roles between the bystander and diquark-participant quarks. 



The Faddeev equation is detailed in Appendix A: Faddeev Equation. With all 
its elements specified, as described therein, the equation can be solved to obtain 
the nucleon's mass and amplitude. Owing to Eq. (IA.36|) . in this calculation the 
masses of the scalar and axial- vector diquarks are the only variable parameters. 
The axial- vector mass is chosen so as to obtain a desired mass for the Z\|j and 
the scalar mass is subsequently set by requiring a particular nucleon mass. 

We have written here of desired rather than experimental mass values because 
it is known that the masses of the nucleon and A are materially reduced by 
pseudoscalar meson cloud effects. This is discussed in detail in Refs. [22l [23] . 
Hence, a baryon represented by the Faddeev equation described above must 
possess a mass that is inflated with respect to experiment so as to allow for an 
additional attractive contribution from the pseudoscalar cloud. As in previous 
work [31 [21] and reported in Table [H we require M^ = 1.18 GeV and m/^ = 
1.33 GeV. The results and conclusions of our study are essentially unchanged 
should even larger masses and a smaller splitting M/\ — Mj^j be more realistic, a 
possibility suggested by Ref. [25] . 



^This is natural because the spin- and isospin-3/2 A contains only an axial-vector diquark. The 
relevant Faddeev equation is not different in principle to that for the nucleon. It is described 
e.g. in Ref. [3]. 



4 Variation of nucleon magnetic moments and radii 

Table 1. Mass-scale parameters (in GeV) for the scalar and axial- vector diquark correlations, 
fixed by fitting nucleon and A masses offset to allow for "pion cloud" contributions [22]. We 
also list Ljjp = —Trmjp, which is the width-parameter in the {qq)jp Bethe-Salpeter amplitude, 
Eqs. (|A.30|) & ()A.31|) : its inverse is an indication of the diquark's matter radius. Row 3 illustrates 
effects of omitting the axial-vector diquark correlation: the A cannot be formed and Mn is 
significantly increased. Evidently, the axial- vector diquark provides significant attraction in the 
Faddeev equation's kernel. 

Mtv Ma mQ+ mi+ ujq+ uji+ 

1.18 LSS 0.796 0.893 0.56=1/(0.35 fm) 0.63=1/(0.31 fm) 

1.46 0.796 0.56=1/(0.35 fm) 



3 Magnetic Moments and Charge Radii 

3.1 Background 

The nucleoli's electromagnetic current is 

MP',P) = ieu{P')A^{q,P)u{P), (3) 

= ieu{P') (^^FiiQ^) + -}-a^,Q,F2{Q^)\ u{P) , (4) 

where P {P') is the momentum of the incoming (outgoing) nucleon, Q = P' — P, 
and Fi and F2 are, respectively, the Dirac and Pauli form factors. They are the 
primary calculated quantities, from which one obtains the nucleon's electric and 
magnetic (Sachs) form factors 

Ge{Q') = Fi(Q') - ^^2(Q') , GMiQ') = F^iQ'') + F2(Q2) . (5) 
The nucleons' magnetic moments are defined through 

lJin = Kn = GIM , /ip = 1 + Kp = G^/O) , (6) 

where kat, N = n,p, are referred to as the anomalous magnetic moments. Of 
course, the nucleon's electric charges, G^(0), are conserved and cannot depend 
on the current-quark mass. That is not true of their electric and magnetic radii: 



r? := -e±Glis) 



^2 
' ' n 



s=0 



-.^Glis) 



K? := -6^^lnGUs) 



(7) 

s=0 

(8) 



=0 

In order to calculate the magnetic moments and charge radii, and their depen- 
dence on current-quark mass, one must know the manner in which the nucleon 
described in Sect. [2] couples to a photon. That is derived in Ref. [26], illustrated 



in Fig. [2] and detailed in Appendix C: Nucleon- Photon Vertex. Naturally, as ap- 
parent in that Appendix, the current depends on the electromagnetic properties 
of the diquark correlations. 
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axial vector scalar 




Figure 2. Vertex which ensures a conserved current for on-shell nucleons described by the 
Faddeev amphtudes, ^ij, described in Sect .[2] and [Appendix A:] Faddeev Equation. The single 
hue represents S(p), the dressed-quark propagator, Sec. lA.2."n and the double line, the diquark 
propagator, Sec. lA.2.^ F is the diquark Bethe-Salpeter amplitude. Sec. I A. 2.^ and the remaining 
vertices are described in |Appendix C:| the top-left image is Diagram 1; the top-right. Diagram 2; 
and so on, with the bottom-right image. Diagram 6. 



3.2 Pseudoscalar meson loops 

The framework we have described hitherto provides what might be caUed the 
quark-core contribution to the nucleons' electromagnetic form factors. As with 
the mass, the nucleons' magnetic moments, and charge and magnetic radii receive 
material contributions from the pseudoscalar meson cloud |27l |28] . There are two 
types of contribution: regularisation-scheme-dependent terms, which are analytic 
functions of m in the neighbourhood of vanishing current-quark mass, m = 0; 
and nonanalytic scheme- independent terms. 

For the magnetic moments and radii the leading-order scheme-independent 
contributions are 1291 



/ \1— loop 
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(9) 
(10) 

(11) 



where, experimentally, gA = 1-26, /^ = 0.0924 GeV= 1/(2.13 fm) and /iy = 
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IJ'P — fJ-n = 4.7. These terms reduce the magnitude of both neutron and proton 
magnetic moments, and increase the magnitudes of the radii. 

While these scheme-independent terms are important, at physical values 
of the pseudoscalar meson masses they do not usually provide the dominant 
contribution to observables. That is provided by the regularisation-parameter- 
dependent terms, as is apparent for bar yon masses in Ref. [22] and for the pion 
charge radius in Ref. [30]. This is particularly significant for the magnetic mo- 
ments, in which connection the regularisation-scheme-dependent terms provide 
a nonzero contribution in the chiral limit and have the net effect of increasing 

This last fact was overlooked in Ref. [24] so that Eq. (82) therein is a poor 
model for the net contribution to fin,p from pseudoscalar meson loops. A minimal 
improvement is 

(/../,)-'°-^ = (<?. ± t^ -^) ^ arctan(^) , (12) 

where /x^^ are the chiral limit values of the meson loop contributions. Equa- 
tion (J12p reproduces Eq. ([9]) but also enables us to express reasonable estimates 
of the regularisation-parameter-dependent parts of the chiral loopsQThe param- 
eter A is a mass-scale. Its presence echoes a physical effect; namely, that finite 
meson size guarantees an intrinsic regularisation of loop integrals. Emulating 
Ref. [M], we employ A = 0.3GeV= l/[0.66fm]. 

A recent estimate from numerical simulations of lattice-regularised QCD [31] 
gives the following chiral-loop contributions to the nucleons' magnetic moments 
at the physical pion mass: 



< = -0.40, /i^ = 0.24. (13) 



They are obtained witqj 



//^ = -1-05, /i^' = 0.88, (14) 



in Eq. (J12p . Equation (J13p suggests that the dressed-quark core we have described 
above should yield the follow magnetic moments for the nucleons at the physical 
current-quark mass, Eq. ()A.23p : 

/,^te9) = -1.51, /i^('?'?) = 2.55. (15) 

These values provide the comparison column in Table [2j 

It was argued in Ref. [2l] that the dressed-quark core described herein is 
compatible with augmentation by a sensibly regulated pseudoscalar meson cloud. 
Indeed, it was argued strongly that agreement with experimental data should not 
be achieved in the absence of such contributions. Given Eq. (J12p . we reconsider 
this below. 



^NB. The expression in Eq. (|12|l vanishes when the pion mass is much larger than the regu- 
larisation scale. This is required because very massive states must decouple from low-energy 
phenomena. 

^Variations in the regularisation mass-scale have no impact; e.g., using A = 0.5 GeV in Eq. (|12|l 
alters these results by less-than 2%. 
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Table 2. Values of the magnetic moments defined in Eq. ([6]) calculated with the diquark mass 



values in Table. [T] and Xi+ = 1- Experimental values 


are |32] /i„ = 


-1.91, /ip = 2.79. 




Eq. (fT^ - quark-core 


Hi+ =0 


^^+ = 0.37 = 


/i+ 


target value 


KT = 


Kr = 0.12 = 


'. Kn~ 


fiT'' -1.51 


-1.44 


-1.51 




M?"^") 2.55 


2.39 


2.55 





3.3 Magnetic properties of axial-vector diquarks 



As explained in Appendix C: Nucleon-Photon Vertex, the nucleons' electromag- 
netic current involves three parameters, which characterise the axial-vector di- 
quarks' electromagnetic properties; viz., the diquarks' magnetic moment - fii+, 
their quadrupole moment ~ Xi+ > and the scalar <-> axial- vector transition strength 
- KT- Calculations have shown [23] that the quadrupole moment has no mate- 
rial impact. Its greatest effect is on the neutron's charge radius, Eq. ([7]). How- 
ever, that is notoriously sensitive to details of the neutron's Faddeev amplitude, 
which determine the magnitude and sign of the Dirac form factor's slope, and to 
meson-related contributions, so that, within the accuracy of the model described 
herein, the net result for all quantities considered is essentially independent of 
Xi+ £ [0, 2]. Variations in the diquarks' magnetic moment and transition strength 
modify the nucleons' radii by < 2%, which is insignificant. Their only real impact 
is on the nucleons' magnetic moments. 

These observations mean that in order to reconsider the claim reiterated 
above, it is only necessary to reanalyse the magnetic moment calculations in 
Ref. [23] . In this case Column 3 of Table [2] reports that with magnetically in- 
ert diquarks the quark-core magnetic moments are underestimated. Moreover, 
Column 4 indicates that there are values of the diquarks' magnetic moment 
and transition strength for which the required quark-core moments, Eq. (J15p . 
can accurately be reproduced. The magnitudes of /ii+ and kt that effect this 
are much smaller than those for an on-shell axial- vector, which were the focus 
of Ref. [231 . This marked suppression is physically reasonable because the di- 
quarks are far from being on-shell within the nucleon, a fact observed explicitly 
in Ref. [33]. Hence, the dominant contributions to the form factors are provided 
by Diagrams 1 and 3. 

It follows that the dressed-quark core defined herein is compatible with dress- 
ing by a sensibly regulated pseudoscalar meson cloud so long as the magnetic 
interactions of axial-vector diquark correlations within the nucleon are commen- 
surate in magnitude with those of the meson cloud. (See also Table O) 

4 Results and Discussion 

We can now describe the method by which we computed the m-dependence of 
the nucleons' magnetic moments and radii. 

We first solved the Faddeev equation at each of a number of current-quark 
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masses within a domain that includes the physical current-quark mass, m^'^^^ 
given in Eq. ()A.23p . In our Faddeev equation the dependence on current-quark 
mass appears explicitly in the dressed-quark propagator described in Sect. lA.2!T] 
and implicitly in the diquark masses. At each value of m the dressed-quark 
propagator is obtained in an obvious manner, Eqs. ()A.19p and (1A.20|) . while the 
diquark masses are given by 

mjP (m) = mfj' + Uqq (m/mP^y*^ - 1) , (16) 

with mF,p^ given in Row 2 of Table [T] and aqq ~ 26MeV being the diquarks' 
cr-term [3l U]. These observations prescribe the m-dependence of the Faddeev 
kernel. Solving the Faddeev equation then provides a range of nucleon Faddeev 
amplitudes, one for each selected current-quark mass. We worked with 

m = 3 MeV + j 6m , 5m = 1 MeV , j = 0, 1, . . . , 7. (17) 

The next step was to calculate, at each value of m in Eq. (J17p . the nucleons' 

electric and magnetic form factors on Q^ € (0,1] GeV^. In doing this we also 

allowed the diquarks' magnetic properties to evolve with current-quark mass. 

jp jp J- jp 

Owing to the fact that F^ := F C', where F is a diquark's Bethe-Salpeter 

amplitude, satisfies approximately the same Bethe-Salpeter equation as a J~^ 

colour-singlet meson )16j . we referred to Ref. [21] for guidance in constraining 

the variation of the diquarks' magnetic parameters. For the p-meson we inferred 

that in the neighbourhood of the physical u/d-quark mass 

/^p("^) = f^f'''" + 0.002 (m/mphys - 1) , /uf ^' = 2.01 , (18) 

from which one extracts the renormalisation group invariant ratio 

5/i„ ,6m 

-^ = 0.001. 19 

/ip m 

Plainly, the response is weak and (perhaps surprisingly, given Table H] herein) the 
/9- meson magnetic moment increases with increasing current-quark mass. Based 
on this analysis, in our calculations we employed 

^^+ (m) = jj.{+ + q'' 0.002 (m/mP^y*^ - 1) , (20) 

Kr{m) = K^ -F ?'' 0.002 (m/mP^y' - 1) , (21) 

where fi^+ and k^- are given in Table [2l Since the response in Eq. (jlSp is weak 
and the sign unexpected, we included the parameters i;^'^ = ±1 and <;^'^ = ±1 in 
Eqs. (j20p . ()2ip so that we could estimate the error associated with employing 
these formulae. 

We judge it worthwhile to dwell a little longer on properties of mesons that 
are relevant to our discourse. From Fig. 7 of Ref. [21] we infer 

6r^ 6m _ 6{r'^f 6m 



/— «^n^/— «-0-06. (22) 



V. V. Flambaum, et al. 9 

Table 3. Quark-core and pseudoscalar loop [Eqs. (|12p . (|26p . (|27p ] contributions to the moments 
and radii, calculated at the physical current-quark mass, Eq. (|A.23|) . The radii are listed in fm^. 
Experimental values are quoted from Ref. [32], where available, and otherwise from Ref. [37| . 





^^n 


I^P 


(r'n) 


{rl) 


{(r^?) 


{(r^?) 


experiment 


-1.91 


2.79 


-(0.34)2 


(0.88)2 


(0.89)2 


(0.84)2 


qigq) 

vr-loop 


-1.51 
-0.40 


2.55 
0.24 


-(0.13)2 
-(0.47)2 


(0.57)2 
(0.47)2 


(0.51)2 
(0.61)2 


(0.50)2 
(0.61)2 


total 


-1.91 


2.79 


-(0.49)2 


(0.74)2 


(0.79)2 


(0.79)2 



This dressed-quark core value possesses the same sign and a similar magnitude 
to the response of the proton's dressed-quark core radii, described subsequently 
in connection with Table [H We have also analysed the response of the pion's 
dressed-quark core. A quick estimate can be obtained via the Nambu-Jona- 
Lasinio model, in which [33] 

It is straightforward to solve this model's gap equation, and therefrom calculate 
/tt and its dependence on current-quark mass. We find thereby 

.^f:=^/^ = -0.05. (24) 



The same result is found by evaluating Eq. (p3|) with the ?7i-dependence of /tt 
depicted in Fig. 3 of Ref. [35], which is obtained from a rainbow truncation of 
QCD's gap equation. Finally, from Fig. 6 of Ref. [36], which reports an internally 
consistent Dyson-Schwinger equation calculation, we infer q^^ = —0.06. A natural 
scale associated with the response of dressed-quark core radii is now apparent. 

Given Eq. (J18p it was not surprising for us to discover that the nucleons' 
magnetic moments respond slowly to changing m. Therefore very precise results 
for the form factors are required in order to compute the slope. In order to 
achieve the numerical accuracy necessary we made significant modifications to 
the computer code employed in Ref. [21] . These changes reduced execution times 
by an order of magnitude, and made practical and reasonable the use of 10* 
adaptive Monte-Carlo sample points in the evaluation of Diagrams 3, 5 and 6, 
which are two- loop integrals. 

Proceeding as described above we arrived at a set of form factors, one at each 
value of the current-quark mass in Eq. (|17p , and thereby the magnetic moments 
and radii as a function of m. Each quantity, designated generically by Q, can be 
expressed in the form 

Q = Qliqq) + g- (25) 

where Q'^^ii' is provided by our Faddeev equation results and Q" is an estimate of 
the contribution from pseudoscalar meson loops Q At the physical current-quark 



*NB. By construction the Faddeev equation is intended to describe a nucleon's dressed-quark 
core. It explicitly excludes all diagrams that can appear in the calculation of Q^ . 
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Table 4. Calculated variation of the nucleons' magnetic moments and radii. We list qq := 
-^/^' where s = q{qq), ^, t = total= q{qq) + n, and fa = -iS^/^. 



0.016 


0.828 


0.186 


JfJ.n 

0.012 


0.029 


1.311 


4 
0.139 


J^lp 
0.026 


y^2 
0.551 


0.477 


0.483 


0.041 


0.014 


0.477 


0.202 


0.008 


0.059 


0.554 


0.350 


0.024 


0.061 


0.554 


0.356 


0.024 



mass, Eq. ()A.23p . our calculated results for the former are given in Table [3j To 
calculate the latter we employed Eq. ([12]) and [Ml [38] 



(^^r^f^i-ioop- = -i±|#ln( ""' ^ 



^^ ' 32^2 j^ ^m2+A2^ 

+T^T2 arctan( ) . (27) 

These values, too, are reported in the Table, from which it is apparent that 
pseudoscalar meson loops contribute materially to the nucleons' static properties. 
As already remarked, the neutron's electric charge radius is particularly sensitive 
to details of the neutron's Faddeev amplitude and, we see here, to the size of the 
meson cloud contribution. In the computation of form factors, that contribution 
remains significant for Q^ < 2GeV2. (See, e.g., Ref. [39].) 

The dependence of the magnetic moments and radii on current-quark mass 
is expressed through 

,e=f^'"" + ^-V™. (28) 

\ dm am I 

We have described in detail the manner in which we calculated the first term. The 
second is computed from Eqs. (|12p . (j26p . (j27p using Eq. ()A.27p . Our results are 
listed in Table l4)i from which it is immediately apparent that the moments and 
radii all decrease with increasing current-quark mass. Moreover, for nucleons, 
which are composed of u/d valence quarks, the response of the contributions 
from the dressed-quark core to a change in current-quark mass is 8 it 3% of that 
arising from the pseudoscalar meson cloud. 

In connection with the magnetic moments it is natural to make a comparison 
with Ref. [12], which reports in Eq. (A7): jIq = —2.58, jIq = 3.48. These values 
can be compared with our results 

t^l^'"'' + /in° = -2.56 , /i^^^-?) + ^xf = 3.43 . (29) 



^NB. The (?^,<;'*)-variation discussed in connection with Eqs. (|20|) . (|2H) affects no result by 
more than 2%. 
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In addition, Eqs. (25), (27) list q^^ = 0.118, q^^ = 0.087 and Eqs. (Al), (A5), 
(A6) enable one to calculate /^„ = 0.010, /^^ = 0.012. Our estimation of the 
loop contributions is cruder than that in Ref. [12j but, nevertheless, our values 
are larger by only a factor of ~ 1.5. Our dressed-quark core contribution to the 
variation of /i^ is almost the same, whereas for the proton it is larger by a factor 
of twoO We emphasise that our results were obtained through an internally 
consistent calculation performed directly at the physical light-quark current- 
mass. In contrast, Ref. [12] extrapolated lattice-regularised QCD results obtained 
at approximately the strange-quark mass. 

These comments notwithstanding, our calculations plainly provide an in- 
dependent confirmation of the magnitude and sign of the effects discussed in 
Ref. [12]. This agreement supports a view that all the results reported in Tabled] 
are reliable estimates. 

5 Epilogue 

We used a Poincare covariant Faddeev equation model for the dressed-quark 
core of the nucleon, augmented by a nucleon-photon vertex which automatically 
fulfills the Ward-Takahashi identity for on-shell nucleons and a rudimentary es- 
timate of the contribution from pseudoscalar meson loops, to obtain insight into 
the response of the nucleons' static electromagnetic properties to changes in 
current-quark mass. 

Our key results are discussed in connection with Table [U and summarised by 
the renormalisation group invariant ratios presented here: 



Q 



O I m 



fJ-n fJ-p rl rl {rf^f (r 



£_ 



0.186 0.139 0.483 0.202 0.350 0.356 



(30) 



Q 

Those for the magnetic moments can assist in constraining the allowed tempo- 
ral variation of the current-quark mass via, e.g., experiments with atomic clocks 
[21 [H] and various astrophysical measurements [2] [13l HO] [H] . Our results also 
suggest that observables dependent on the nucleons' magnetic and charge radii 
might provide a useful means by which to place limits on the allowed variation in 
Nature's fundamental parameters. It is notable, for example, that the calculated 
energy levels and transition frequencies in Hydrogen and Deuterium, which are 
some of the most precise theoretical predictions in physics, are quite sensitive to 
the value of the proton's charge radius [52]. Moreover, all hyperfine transition fre- 
quencies in atomic clocks [21 [l3] and astrophysics [21 [131 BOl W\\ react to a change 
in nucleonic charge and magnetisation distributions: the former alters electronic 
wave functions and the latter changes the hyperfine interaction Hamiltonian. 

As a byproduct of this study, we arrived at an improved understanding of 
diquark correlations. Within the nucleon they are usually far from being on- 
shell and hence it is a poor approximation to represent, for example, the active 
magnetic properties of the axial- vector correlations by on-shell values. 
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Appendix A: Faddeev Equation 

A.l General structure 

The nucleon is represented by a Faddeev amplitude 

\lr ^ tpi + >I'2 + 'Pa , (A.l) 

where the subscript identifies the bystander quark and, e.g., !f'i,2 are obtained from •f's by a 
correlated, cyclic permutation of all the quark labels. We employ the simplest realistic repre- 
sentation of 'P. The spin- and isospin-1/2 nucleon is a sum of scalar and axial-vector diquark 
correlations: 

1'3{p,,m,T,)=J^t +J^t, (A.2) 

with {pi,ai,Ti) the momentum, spin and isospin labels of the quarks constituting the bound 
state, and P — pi + P2 + P3 the system's total momentum[j 
The scalar diquark piece in Eq. l|A.2p is 

f^^ {p,,m, n) = [r°* (ip[i2] ; K)]l\X ^°^ (A") [5(£; P)n{P)]ll , (A.3) 

where: the spinor satisfies ( [Appendix Ji:\ Euclidean Conventions) 

{i'yP + M) u{P) = = u{P) (i-y- P + M) , (A.4) 

with M the mass obtained by solving the Faddeev equation, and it is also a spinor in isospin 
space with ip+ = col(l,0) for the proton and ip- = col(0, 1) for the neutron; K = pi + p2 ~' 
P{i2}, P[i2] = Pi — P2, a := (— P{i2} + 2p3)/3; Zi" is a pseudoparticle propagator for the scalar 
diquark formed from quarks 1 and 2, and P'^ is a Bethe-Salpeter-like amplitude describing their 
relative momentum correlation; and 5, a 4x4 Dirac matrix, describes the relative quark-diquark 
momentum correlation. (5, P" and Zi" are discussed in Sect. lA. 21 ') The colour antisymmetry 
of 1^3 is implicit in P , with the Levi-Civita tensor, ecic2C3, expressed via the antisymmetric 
Gell-Mann matrices; viz., defining 

{H^ =i\\H'^ = ~i\\H'^ = i\^}, (A.5) 

then ecic2C3 = {H'"')cic2- [See Eqs. COOJI . (|A3T|| .] 
The axial- vector component in Eq. (IA.2|I is 

M'^ip,,a,,n) = lt'P^^{lpi^2y,K)]l\XAlt{K)[Al{£;P)u{P)]ll, (A.6) 

where the symmetric isospin-triplet matrices are 

t+ = ^(r° + r^),t"=r\t- = -L(r°-r^), (A.7) 



and the other elements in Eq. (|A.6|) are straightforward generafisations of those in Eq. HA.3[I . 

The general forms of the matrices S{£;P), A]^{1; P) and 'D^p{£; P), which describe the 
momentum space correlation between the quark and diquark in the nucleon are described in 
Refs. [431 144) . The requirement that S{£;P) represent a positive energy nucleon; namely, that 
it be an eigenfunction of A+{P), Eq. (|B.7|) . entails 



S{£;P)^si{£;P)Io+(i-f£-£-PIu) S2{£;P), (A.8) 



^NB. Herein we assume isospin symmetry of the strong interaction; i.e., the u- and d-quarks 
are indistinguishable but for their electric charge. This simplifies the form of the Faddeev 
amplitudes. 
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where (I-D)rs = Srs, i = 1, P = —1. In the nucleon rest frame, si,2 describe, respectively, the 
upper, lower component of the bound-state nucleon's spinor. Placing the same constraint on 
the axial- vector component, one has 



Al{e; P) = J2 P^(^' ^) ^5 A^il; P) , i = +, 0, 



(A.9) 



where 



-f^-PP,,7^ = 7, -f7-PP, 



Ai=7- 
At =i£ 



I' 1 



Al = i7.^7 • 



■At 



One can now write the Faddeev equation satisfied by ^z as 



S{k-P)u{P) 
Al{k-P)u{P) 



(2^) 



■M{kJ-P) 



S{t,P)u[P) 
AU£;P)u{P) 



The kernel in Eq. (|A.ll|l is 



Mik,e;P) 



Moo 
(Mio)p 



{Moi)i 



with 



Moo = P" (fcg - ^g,/2; 4,) 5'^(^,, - fc,) P° {£, - fc„/2; -A:,,) ^(f,) A" {£,,) , 



(A.IO) 



(A.ll) 



(A.12) 



(A.13) 



where: Iq =l + P/3, kg = fc-f P/3, £qq = -£. + 2P/3, kqq = -fc + 2P/3 and the superscript "T" 
denotes matrix transpose; and 



XO y^qq — kq) 1 [Zq — Kqq / Z\ —f^qq) J\^q) ^^v \^qq) -, 

{Mioy^ = r°\kq-£qj2-lqq) 

Xo ("C-gg— /Cgjt 1^^ [£q — Kqq / Z] — Kqq ) J [iq ) Z\ (-^ggj; 

{Mii)% = t'r^\kq ^ lqq/2;£qq) 

X J \^qq ^q) "^ -'- fi \^q <^qql Z\ rZqq) '^{yqj ^pv {^qq) • 



(A.14) 
(A.15) 
(A.16) 



A. 2 Kernel of the Faddeev equation 

To complete the Faddeev equations, Eq. ()A.1H) . one must specify the dressed-quark propagator, 
the diquark Bethe-Salpeter amplitudes and the diquark propagators. 



A. 2.1 Dressed-quark propagator 

The dressed-quark propagator has the general form 



Sip) 



-i-y ■ pav{p ) + (ys{p ) = 1/[J7 ■ pMp ) + B{p 



(A.17) 



and can be obtained from QCD's gap equation; namely, the Dyson-Schwinger equation for the 
dressed-quark self-energy [45]. The gap equation has been much studied and features of its 
solution elucidated. Indeed, it is a longstanding prediction of DSE studies in QCD that for 
light-quarks the wave function renormalisation and dressed-quark mass: 

Z{p^) = 1/A{p^) , Mip^) = Bip^)/A{p^) , (A.18) 

respectively, receive strong momentum-dependent corrections at infrared momenta [451 1461 147| : 
Z{p^) is suppressed and M(p^) enhanced. These features are an expression of dynamical chiral 
symmetry breaking (DCSB) and, plausibly, of confinement [48]. The enhancement of M{p^) is 
central to the appearance of a constituent-quark mass-scale and an existential prerequisite for 
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Goldstone modes. The mass function evolves with increasing p^ to reproduce the asymptotic 
behaviour familiar from perturbative analyses, and that behaviour is unambiguously evident 
for V^ ^ 10 GeV^ [49]. These DSE predictions are confirmed in numerical simulations of lattice- 
regularised QCD [50] , and the conditions have been explored under which pointwise agreement 
between DSE results and lattice simulations may be obtained [511 1521 [55] . 

The impact of this infrared dressing on hadron phenomena has long been emphasised [54] 
and, while numerical solutions of the quark DSE are now readily obtained, the utility of an 
algebraic form for S(p) when calculations require the evaluation of numerous multidimensional 
integrals is self-evident. An efficacious parametrisation of S(p), which exhibits the features 
described above, has been used extensively in hadron studies [55]. It is expressed via 

as{x) = 2rhr{2{x + m))+F{hix)J^{h3x)[ho + h2J^{f.x)\ , (A. 19) 

av{x) = ,^ „ [1 - :F{2ix + m))] , (A.20) 

X + m 



with X = p /X , m — m/X, 



X 



T{x) = ^^^ , (A.21) 



o"s(a;) = Xas{p ) and av{x) = X av{p )■ The mass-scale, X = 0.566 GeV, and parameter 
valueflfl 

fh bo 6i 62 &3 



0.00897 0.131 2.90 0.603 0.185 ' 



(A.22) 



were fixed in a least-squares fit to light-meson observables [561 157] . The dimensionless u = d 
current-quark mass in Eq. ([A.22|) corresponds to 

m = 5.08 MeV =: m^^^'' . (A.23) 

The parametrisation yields a Euclidean constituent-quark mass 

Mi^d = 0.33 GeV, (A.24) 

defined as the solution of p^ — M^{p^). 

The ratio M^ /m = 65 is one expression of DCSB in the parametrisation of S{p). It 
emphasises the dramatic enhancement of the dressed-quark mass function at infrared momenta. 
Another is the chiral-limit vacuum quark condensate [53] 

3 bo ,„ C' 



('^'^)? = A^,^^ln^ (A.25) 



47r2 6163 '" A"^ 



which assumes the value (tIqcd = 0.2 GeV) 

- (w)?=i GeV = (0.221 GeV)^ (A.26) 

A detailed discussion of the vacuum quark condensate in QCD can be found in Ref. [581 [55] 

An exact formula for pseudoscalar meson masses was derived in Ref. [BD] and in the present 
context it can be expressed for the pion as 

firnl = -2m(gg)7GcV , (A. 27) 

where (qq) iccv is an in- pion condensate [61]. In order to calculate this quantity and /^ one 
needs the pion's Bethe-Salpeter amplitude. An Ansatz compatible with the parametrisation of 
the dressed-quark propagator described above is explained in Ref. [57], and together they yield 

- (ot)"gov = (0.250 GeV)^ , ,U = 0.090 GeV , m^ = 0.140 GeV . (A.28) 



*The parameters 60,1,2,3 are assumed to be m-independent. In the current application, that is 
possibly a weakness of the parametrisation. For example, it leads to a constituent-quark cr-term 
that is 30% smaller than that obtained from the solution of a well-constrained gap equation 
[3]. e = 10"'* in Eq. (|A.19[) acts only to decouple the large- and intermediate-p^ domains. 
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A. 2. 2 Diquark Bethe-Salpeter amplitudes 

The rainbow-ladder DSE truncation yields asymptotic diquark states in the strong interaction 
spectrum. Such states are not observed and their appearance is an artefact of the truncation. 
Higher-order terms in the quark-quark scattering kernel, whose analogue in the quark-antiquark 
channel do not much affect the properties of vector and flavour non-singlet pseudoscalar mesons, 
ensure that QCD's quark-quark scattering matrix does not exhibit singularities which corre- 
spond to asymptotic diquark states [62] . Nevertheless, studies with kernels that do not produce 
diquark bound states, do support a physical interpretation of the masses, ™.(gq) p , obtained 
using the rainbow-ladder truncation: the quantity l^ggf ^ — l/jnjgq) ^ rnay be interpreted as a 
range over which the diquark correlation can persist inside a baryon. These observations moti- 
vate an Ansatz for the quark-quark scattering matrix that is employed in deriving the Faddeev 
equation: 

[M,,(fc,Q;if)]^:= J2 r'\k--K)A'\K)r'\q-K). (A.29) 

JP=0+,1+,... 

One means of specifying the F'' in Eq. HA.29|I is to employ the solutions of a rainbow-ladder 
quark-quark Bethe-Salpeter equation (BSE), as e.g. in Refs. [18II19J . Using the properties of the 
Gell-Mann matrices one finds easily that Fq := F^ C^ satisfies exactly the same equation as 
the J~^ colour-singlet meson hut for a halving of the coupling strength [16]. This makes clear 
that the interaction in the 3c (gg) channel is strong and attractive. 

A solution of the BSE equation requires a simultaneous solution of the quark-DSE. However, 
since we have already chosen to simplify the calculations by parametrising S{p), we also employ 
that expedient with F , using the following one-parameter forms: 

F°^{k;K) = _i_J/"C'n5ir2J^(fcV^o+)> (A.30) 

t'F^^{k;K) = --i^iJ%^,,Ct'.F(fcVc^?+), (A.31) 

with the normalisation, A/""^ , fixed by requiring 



o -\ K =-m. 



aQ.^(^'^) 



d^g 



(A.32) 



Q=K 



-,T/ 



n(K,Q) = trJj^F{q-~K)S{q + Q/2)F(q-K)S'{-q + Q/2). (A.33) 

The Ansatze of Eqs. (|A.30|l . HA.3H) retain only that single Dirac-amplitude which would 
represent a point particle with the given quantum numbers in a local Lagrangian density. They 
are usually the dominant amplitudes in a solution of the rainbow- ladder BSE for the lowest 
mass J^ diquarks [17II18J and mesons [611 1631 [64] . 



A. 2. 3 Diquark propagators 

Solving for the quark-quark scattering matrix using the rainbow-ladder truncation yields free 
particle propagators for A''' in Eq. ([A.29[l . As already noted, however, higher-order contri- 
butions remedy that defect, eliminating asymptotic diquark states from the spectrum. The 
attendant modification of A'' can be modelled efficiently by simple functions that are free- 
particle- like at spacelike momenta but pole- free on the timelike axis [62]; namely|j 

A'>\k) = ^^(7^V^o+), (A.34) 

KUk) = (s,, + ^^) ^:F{K'/coI+), (A.35) 



^These forms satisfy a sufficient condition for confinement because of the associated violation 
of reflection positivity. See Sect. 2 of Ref. [48] for a brief discussion. 
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where the two parameters rrijp are diquark pseudoparticle masses and uijp are widths charac- 
terising r''' . Herein we require additionally that 



dK'2 



■ ^(i^V^jp) 



K2=0 



= 1 ^ Lo'^jP = - m^jP 



(A.36) 



which is a normalisation that accentuates the free-particle-like propagation characteristics of 
the diquarks within the hadron. 



Appendix B: Euclidean Conventions 

In our Euclidean formulation: 

4 



{7m,7^} = 2(5p^ ; 7p = 7m; cr^^ = -^biJ.,!!^]; tr[757M7^7p7a] = -4 6^,^^^ ,£1234 = 1. 
A positive energy spinor satisfies 

u{P, s) (J7 • P -I- M) =0 = {i^- P + M) u{P, s) , 
where s = ± is the spin label. It is normalised: 

u{P, s) u{P, s) = 2M 
and may be expressed explicitly: 



j.{P, s) = VM - i£ (T P 



M -i£' 



with £ = iV-P^ + M^ 



X+ 



For the free-particle spinor, u(P, s) = u{P, 3)^74. 

The spinor can be used to construct a positive energy projection operator: 

s — ± 

A negative energy spinor satisfies 

v{P, s) {i-f-P- M) =0 = {i-f P-M) v{P, s) , 



(B.l) 
(B.2) 

(B.3) 

(B.4) 

(B.5) 
(B.6) 

(B.7) 
(B.8) 



and possesses properties and satisfies constraints obtained via obvious analogy with u{P, s) 
A charge-conjugated Bethe-Salpeter amplitude is obtained via 



r{k-P) =C'' r{-k;PfC, 



(B.9) 



where "T" denotes a transposing of all matrix indices and C = 7274 is the charge conjugation 
matrix, C^ = — C. 



Appendix C: Nucleon-Photon Vertex 

In order to explicate the vertex depicted in Fig. [2] we write the scalar and axial- vector compo- 
nents of the nucleons' Faddeev amplitudes in the form [cf. Eq. HA.llf) ] 



^{k;P) 



W^ik-P) 
%{k;P) 



S{k;P)u{P) 
AUk;PMP) 



i = l, 



(C.l) 
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For explicit calculations, we work in the Breit frame: P^ — P/^^ — Qn/2, P^ = P^^ + Qi_i/2 
and P^^ = (0, 0, 0,i^M^ + Q^/4), and write the electromagnetic current matrix element as 
[cf. Eq. (0] 



p'\j;r\p) = A+{p') 



j4„ j4 



"ItiGs + Mn-fy2 — (Ge — Gm) 



p2 



yl+(P) (C.2) 



, = r?P + fc , 


Pq = 'nP' +p, 


d^ fjP - k, 


Pd = vP' - P , 



= l^,^n-p,p')jrip,p'-,k,pmk,p). (c.s) 

In Fig.[2]we have broken the current, J"^^{p, P'\ k, P), into a sum of six terms, each of which we 
subsequently make precise. NB. Diagrams 1, 2 and 4 are one-loop integrals, which we evaluate 
by Gaufiian quadrature. The remainder. Diagrams 3, 5 and 6, are two-loop integrals, for whose 
evaluation Monte-Carlo methods are employed. 

C.l Diagram 1 

This represents the photon coupling directly to the bystander quark. It is expressed as 

jr = S{p,)r'^-{p,;k,)S{k,)(^A"^{ks)+A'^{ks)){27v)^5\p~k~fiQ), (C.4) 

where r^'^[pq; kq) = Qq r^{pq; kq), with Qq = diag[2/3, —1/3] being the quark electric charge 
matrix, and ri^{pq; kg) is the dressed-quark-photon vertex. In Eq. (|C.4p the momenta are 

(C.5) 

with rj + fi = 1. The results reported herein were obtained with rj = 1/3, which provides a single 
quark with one-third of the baryon's total momentum and is thus a natural choice. Notably, 
as our approach is manifestly Poincare covariant, the precise value is immaterial so long as 
the numerical methods preserve that covariance. Calculations converge most quickly with the 
natural choice. 

It is a necessary condition for current conservation that the quark-photon vertex satisfy 
the Ward-Takahashi identity: 

Q^ir^{e.i,i>2) = ^"'(^i) - s-\e.2) , (c.6) 

where Q = ^i — ^2 is the photon momentum flowing into the vertex. Since the quark is dressed. 
Sec. I A. 2. n the vertex is not bare; i.e., Pp(^i,i'2) 7^ 7^- It can be obtained by solving an inho- 
mogeneous Bethe-Salpeter equation, which was the procedure adopted in the DSE calculation 
that successfully predicted the electromagnetic pion form factor |20l 164) . However, since we 
have parametrised S(p), we follow Ref. [SJ and write [BS] 

ir^{h,l2) = iSA{ti,ti) 7m + 2fc^ [n • k^ Aa^hIaI) + AB{i\,il)] ; (C.7) 

with fc = (^1 -h ^2)/2, Q = (^1 - h) and 

^F(^?,^i) = i[P(^?)+P(£i)], AF{i\,il)='^^^^^^^, (C.S) 

where F = A,B\ viz., the scalar functions in Eq. HA.17|) . It is critical that P^j in Eq. (|C.7|) satisfies 
Eq. (|C.6|) and very useful that it is completely determined by the dressed-quark propagator. 



C.2 Diagram 2 

This figure depicts the photon coupling directly to a diquark correlation. It is expressed as 

4" = A\pd)\f^''{pd;kd)r A'{kd)S{kq){27TfS*{p-k + rjQ) (C.9) 



with [f^^iPd-, kd)Y' = diag[Qo+P°^ , Qi+r^^], where Qo+ = 1/3 and P^^ is given in Eq. (ICT2ll . 
and Qi+ = diag[4/3, 1/3, —2/3] with P^ given in Eq. (|C.14|) . Naturally, the diquark propaga- 
tors match the line to which they are attached. 
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In the case of a scalar correlation, the general form of the diquark-photon vertex is 

r°Uli,£2) = 2k^ f+{k\k ■ Q,Q^) + Q^ f-{k\k ■ Q,Q^) , (CIO) 

and it must satisfy the Ward-Takahashi identity: 

Q^r°\£^,£2) = n'-'\£l)-n'>^i£l), n'" (i^) = {A-'" i£^)}-\ (c.ii) 

The evaluation of scalar diquark elastic electromagnetic form factors in Ref. [19] is a first step to- 
ward calculating this vertex. However, in providing only an on-shell component, it is insufficient 
for our requirements. We therefore adapt Eq. (|C.7|l to this case and write 



rf{£i,£2)^k^A^,+ {£l£l), (C.12) 

with the definition of A q+{£i,£2) apparent from Eg. (jC.8|l . Equation HC.12|| is the minimal 
Ansatz that: satisfies Eq. (|C.ll[l : is completely determined by quantities introduced already; 
and is free of kinematic singularities. It implements /_ = 0, which is a requirement for elastic 
form factors, and guarantees a valid normalisation of electric charge; viz., 

lim^rf {£',£) = 2£,^n''\e) '"=°2£, , (C.13) 

owing to Eq. HA.36|) . NB. We have factored the fractional diquark charge, which therefore ap- 
pears subsequently in our calculations as a simple multiplicative factor. 

For the case in which the struck diquark correlation is axial-vector and the scattering is 
elastic, the vertex assumes the form [66j:r°l 

3 

rlU{£i,£2) ^ -Y.rf^is{£i,£2) , (C.14) 

1=1 

^Nit\i{T^p(£)^S^p-£^£fi/£^) 

r!!U^i,i2) = {,£i+£2)^T^x{£i)Txp{,£2)F^(£l,£l), (C.15) 

^l(^i'^2) = [T^c,{£i)Tpp{^2)£ip + T^p{£2)T^p{£i)£2AF2(£\,£l), (C.16) 



rllpi^i,i2) = -—^i^i+^2)^.Tc,p{£i)£2pTi3x{£2)£ixFa{£i,£i). (C.17) 

This vertex satisfies: 



eic.r'^l0{£i,£2) = o^r'^t0{eij2)£2p, (c.is) 

which is a general requirement of the elastic electromagnetic vertex of axial- vector bound states 
and guarantees that the interaction does not induce a pseudoscalar component in the axial- 
vector correlation. We note that the electric, magnetic and quadrupole form factors of an 
axial- vector bound state are expressed [66] 

2 



Cr (Q') = Fi + i r,^ G'q (Q^) , n+ = -^ (C.19) 



G^(Q') = -i^2(Q') , (C.20) 

Gq"'(0') =Fi(0')+F2(Q')-f(H-ri+)i^3(Q'). (C.21) 

Owing to the fact that Fq '■— F' C' satisfies exactly the same Bethe-Salpeter equation as 
the J~^ colour-singlet meson but for a halving of the coupling strength, the vector meson form 
factor calculation in Ref. [21] might become useful as a guide in understanding the form factors 
in Eqs. (|C.14|l - (|C.17|l . However, in providing only an on-shell component, that information is 
insufficient for our requirements. Hence, we employ the following Ansatze: 

Fi{£l£l) = A^,^{£\,£l), (C.22) 

F2[,£\AI) = -Fi + (l-ri+)(ri+Fi + l-Mi+)rf(^i+) (C.23) 

fM.£%) = -(xi+(l-ri+)d(Ti+)+Fi+F2)d(ri+), (C.24) 



^°If the scattering is inelastic the general form of the vertex involves eight scalar functions [67] 
For simplicity, we ignore the additional structure in this Ansatz. 
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with d(x) — 1/(1 + x)^. This construction ensures a valid electric charge normalisation for the 
axial- vector correlation; viz., 

Urn r'^t^{£',£) = T^p{£) -^ n'* (£^) ' = ° T^p(£)2£^ , (C.25) 

owing to Eq. (|A.36|) . and current conservation 

lim Qf.r^lB{£i,£2)^0. (C.26) 

The diquark's static electromagnetic properties follow: 

Gf (0) = 1 , G'm (0) = Mi+ , G^a (0) = -Xi+ • (C.27) 

For an on-shell or pointlike axial-vector: fii+ = 2; and Xi+ = li which corresponds to an 
oblate charge distribution. In addition, Eqs. (|C.14p - ()C.17|) with Eqs. HC.22f| - (|C.24|l realise the 
constraints of Ref. [68] ; namely, independent of the values of /ii+ & Xi+ 1 the form factors assume 
the ratios 

GI\q') : G'm{Q') ■■ G'q{Q^) '^'=°° (1 - |ri+) : 2 : -1. (C.28) 

We note that within a nucleon the diquark correlation is not on-shell. Hence, in con- 
trast with Ref. [25, herein we do not assume that point-particle values for the magnetic and 
quadrupole moments in Eqs. (|C.27|I serve as a good reference point. For the processes described 
by Fig. [21 the values can be much smaller in magnitude [33) . as we find in Table [51 

C.3 Diagram 3 

This image depicts a photon coupling to the quark that is exchanged as one diquark breaks up 
and another is formed. It is expressed as 

^r = -lsik,)A\k,)r\pi,ka)S^{q)rr^iq',q)S^{q')r'^ip'2,p,)AHp<i)S{p,) , (C.29) 



wherein the vertex r)?" appeared in Eq. (|C.4|I . While this is the first two- loop diagram we have 
described, no new elements appear in its specification: the dressed-quark-photon vertex was 
discussed in Sec. IC.ll In Eq. HC.29|) the momenta are 

q = 77P — rjP' — p — k , q' — rjP' — rjP — p — k , 

Pi = {p,-q)/2, ' p'2 = {^k, + q')/2, (C.30) 

pi = (P, - 9')/2 , p2^{-kg + q)/2. 

It is noteworthy that the process of quark exchange provides the attraction necessary in 
the Faddeev equation to bind the nucleon. It also guarantees that the Faddeev amplitude has 
the correct antisymmetry under the exchange of any two dressed-quarks. This key feature is 
absent in models with elementary (noncomposite) diquarks. The full contribution is obtained 
by summing over the superscripts i,j, which can each take the values 0^, 1^. 

C.4 Diagram 4 

This differs from Diagram 2 in expressing the contribution to the nucleons' form factors owing 
to an electromagnetically induced transition between scalar and axial- vector diquarks. The 
explicit expression is given by Eq. (|C.9|) with [r^'^{pd\ fcd)]'"^ = 0, and [r^'^{pd\ ^d)]^'^ = Psa 
and [r^''{pd\kdy^'^ = Fas- This transition vertex is a rank-2 pseudotensor, kindred to the 
matrix element describing the p7*7r'' transition [69], and can therefore be expressed 

rrA{^u£2) = -rrs{£i,£2) = j^T{£^,£2)e^o,px£ip£2x, (C.31) 

where 7, a are, respectively, the vector indices of the photon and axial-vector diquark. For 
simplicity we proceed under the assumption that 

T{£l,£2)=^^T■, (C.32) 
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viz., a constant, for which a typical on-shell value is kt ~ 2 [70]. However, as with /ii+ and 
Xi+ , we recognise herein that this value is not a useful reference point because, for the processes 
described by Fig.[2j kt can be much smaller in magnitude. 

In the nucleons' rest frame, a conspicuous piece of the Faddeev amplitude that describes 
an axial-vector diquark inside the bound state can be characterised as containing a bystander 
quark whose spin is antiparallel to that of the nucleon, with the axial-vector diquark's spin 
parallel. The interaction pictured in this diagram does not affect the bystander quark but the 
transformation of an axial-vector diquark into a scalar effects a flip of the quark spin within the 
correlation. After this transformation, the spin of the nucleon must be formed by summing the 
spin of the bystander quark, which is still aligned antiparallel to that of the nucleon, and the 
orbital angular momentum between that quark and the scalar diquark. This argument, while 
not sophisticated, does motivate an expectation that Diagram 4 will impact upon the nucleons' 
magnetic form factorsr^l 

C.5 Diagrams 5 & 6 

These two-loop diagrams are the so-called "seagull" terms, which appear as partners to Dia- 
gram 3 and arise because binding in the nucleons' Faddeev equations is effected by the exchange 
of nonpointlike diquark correlations [26]. The explicit expression for their contribution to the 
nucleons' form factors is 






- r\pi,k^)S^{q)Xf,{~k„ -q,Pd)) A^pa)S{p,) , (C.33) 

where, again, the superscripts are summed. 

The new elements in these diagrams are the couplings of a photon to two dressed-quarks as 
they either separate from (Diagram 5) or combine to form (Diagram 6) a diquark correlation. As 
such they are components of the five point Schwinger function which describes the coupling of a 
photon to the quark-quark scattering kernel. This Schwinger function could be calculated, as is 
evident from the computation of analogous Schwinger functions relevant to meson observables 
[71j . However, such a calculation provides valid input only when a uniform truncation of the 
DSEs has been employed to calculate each of the elements described hitherto. We must instead 
employ an algebraic parametrisation [21], which for Diagram 5 reads 

Xf {k, Q) = e^y ^f'Q^Q2 [^'"(fc - Q/2) - r'"{k)] 



Ak-Q-Q2 



^|^[r-^^(/c + Q/2)_r^w], (C.34) 



with k the relative momentum between the quarks in the initial diquark, eby the electric charge 
of the quark which becomes the bystander, and ecx the charge of the quark that is reabsorbed 
into the final diquark. Diagram 6 has 

xfik,Q) = e,, ^f:'^;^' [r-^"(fc + Q/2) - r-^"(fc)] 



-I- eo 



4fc-Q + 
4fc^ - Qf. 



Ak-Q-Q 



^ \r'\k - Q/2) - r'\k)] , (C.35) 



where F {(.) is the charge-conjugated amplitude, Eq. (|B.9[I . Plainly, these terms vanish if the 
diquark correlation is represented by a momentum-independent Bethe-Salpeter-like amplitude; 
i.e., the diquark is pointlike. 



^'^ Another component of the amplitude has the bystander quark's spin parallel to that of the 
nucleon while the axial- vector diquark's is antiparallel: this q^ © (??);+ system has one unit of 
angular momentum. That momentum is absent in the 5^ © (q?)o+ system. Other combinations 
also contribute via Diagram 3 but all mediated processes inevitably require a modification of 
spin and/or angular momentum. An analysis of the contribution from quark orbital angular 
momentum to a nucleon's spin is presented in Ref. [44) . 
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It is naturally possible to use more complicated Ansdtze. However, like Eq. (|C.12[l . 
Eqs. (|C.34[l & (|C.35[l are simple forms, free of kinematic singularities and sufficient to ensure 
the nucleon-photon vertex satisfies the Ward-Takahashi identity when the composite nucleon 
is obtained from the Faddeev equation. 

Appendix D: Charge Symmetry 

Our analysis assumes rriu ~ rud- Hence the only difference between the u- and rf-quarks is 
their electric charge. Our equations and computer codes therefore exhibit the following charge 
symmetry relations: 

where /xJJ means the contribution from the u-quark to the magnetic moment of the neutron, 
etc.; and furthermore 

SfJ-n = -25/ip , Sfit = ~-Sfj.p , (D.2) 

where, in this section, Sfi^ means the variation in ^^ owing to a small change in current-quark 
mass. Using these equations, one obtains 

Sfip _ 5^p + S^p 5 fin Sfip + 4(5 jUp 

Our Faddeev equation model yields 

Mp=2.40, <5^;j = -0.032 , 

fip = 0.15, Sfip = 0.0034, ^ ■ ' 

from which we obtain the results in Tables |3] and ID 

It is interesting to provide a context for the results in Eq. (|D.4[l . Suppose one were required 
to reproduce Hp in Eq. (|15|l with nonrelativistic pointlike constituent-quarks. Such quarks 
have the magnetic moments: 

^J'U = '2fJ.Q , ^^D = -ifiQ , (D-5) 

in terms of which 

QiQQ) 4 1 u , D /Ti fl^ 

Mp = gMc - gMo =: tJ-p + fJ-p ■ (D-6j 

With fiQ =0.85 one reproduces Eq. (|15|) and finds 

Up =2.27, Up =0.28. (D.7) 



A comparison between Eqs. (|D.4[I and (|D.7|l indicates the presence of correlations in our Faddeev 
amplitude for the nucleon. Relative to a generic constituent-quark model, they increase the 
probability for a ii-quark to have its spin aligned with that of the proton, and markedly decrease 
that probability for the rf-quark. 
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